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CN ■ ABSTRACT 

We revisit the shapes of isophotes for elliptical (E) and lenticular (SO) galaxies by 
studying 847 nearby early-type galaxies selected from the Sloan Digital Sky Survey 
Data Release 4 with velocity dispersions above 200 kms" 1 . The IRAF task ellipse 
was used to derive the deviations of the isophotes from pure ellipses (Fourier coeffi- 
cients a^/a and a^/a), position angles and ellipticities as a function of radius. We show 
' the statistical distributions of the 0,3/0, and 04/0 parameters as a function of velocity 

, dispersion, magnitude, and colour. The a 4 /a parameter is correlated with the elliptic- 

ity and absolute Petrosian magnitude of galaxies. No significant correlation was found 
between the 04/ a parameter with colour and velocity dispersion. A cross-correlation 
between the FIRST survey and the SDSS data reveals a lack of powerful radio emitters 
in discy E/SOs, as previously found by Bender et al. We also find that boxy E/SOs 
\ favor denser environments while discy E/SOs favor more isolated environments. The 

median values of changes in the ellipticity and position angle between one and one 
and a half Petrosian half-light radii in the isophotes are about —0.023 and 1.61° re- 
spectively. The average change in the position angle is much larger, about 4.12°, due 
to an extended tail. The change in ellipticity is weakly correlated with the ellipticity 
' itself, with an increasing ellipticity for galaxies with higher ellipticity as the radius 

increases. The isophote parameters for the 847 galaxies are available online. 
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1 INTRODUCTION 

The morphologies of galaxies carry important information 
on how they form, evolve and how they are supported dy- 
namically. For example, disc galaxies are highly flattened 
systems that are rotationally supported. It implies that an- 
gular momentum must have played a crucial role in de- 
termining their properties, such as sizes and the Tully- 
Fisher relation (e.g., Mo, Mao & White 1998 and references 
therein). Similarly, shapes of elliptical (E) and lenticular 
(SO) galaxies carry important information about their for- 
mation history and dynamical state. 

The isophotes of elliptical galaxies and SO's can be well 
htted by ellipses to the zeroth order. Small but significant de- 
viations from perfect ellipses exist. Around 1980's, CCDs be- 
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came more widely used in astronomy and they enabled accu- 
rate photometric measurements of galaxies. Several authors 
studied the isophotal shapes of elliptical galaxies (Lauer 
1985; Bender & Mollernhoff 1987; Jedrzejewski 1987). In an 
influential study, Bender et al. (1988) presented the isopho- 
tal analysis of a sample of 69 nearby bright elliptical galax- 
ies. This study used CCD photometry in the Cousins V, 
R, I under typical seeing conditions of about 2 arcseconds 
with a field of view of 3 by 4 arc minutes. Combining CCD 
data from other studies, Bender et al. (1989) classified 109 
early-type galaxies according to the isophotal shapes using 
Fourier expansions in the polar angle (see ^2.2.21 . It was 
found that the most significant non-zero component of the 
Fourier analysis is the 0,4 parameter (corresponding to the 
cos 40 term). Using the sign of this parameter, E/S0 galax- 
ies were classified into discy (04 > 0) and boxy (04 < 0) 
galaxies. Bender et al. (1989) showed that this parameter 
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has significant correlations with the radio and X-ray prop- 
erties in elliptical galaxies. It is now increasingly clear that 
early- type galaxies can be subdivided into two classes: discy 
ellipticals tend to be fainter, rotationally supported, lack X- 
ray and radio activities and have power-law nuclear light 
profiles while boxy ellipticals tend to be brighter, supported 
by random motions, have significant X-ray and radio activi- 
ties and cored nuclear profiles (Ferrarese et al. 1994; van den 
Bosch et al. 1994; Lauer et al. 1995; Faber et al. 1997; Rest et 
al. 2001; Lauer et al. 2005). In addition, the shapes of ellip- 
tical galaxies may also be correlated with their ages (Ryden, 
Forbes & Terlevich 2001). Pellegrini (1999, 2005) updated 
the analysis of the correlation between galaxy shapes and 
X-ray properties for early-type galaxies using ROSAT and 
Chandra data. Pasquali et al. (2006) found, among other 
things, that the percentages of discy and boxy ellipticals at 
z ~ 1 are similar to the local values, and their characteristic 
shapes also follow the same correlations of local ellipticals, 
with an exception of the 0,3/0, parameter which appears to 
be larger than that for the local counterparts. 

The ellipticities and position angles in ellipticals vary 
as a function of distance from the galaxy centre. Di Tullio 
(1978, 1979) found that for isolated ellipticals, on average, 
the ellipticity decreases as a function of radius, while for 
ellipticals in clusters and groups, the ellipticity either in- 
creases, decreases or peaks as a function of radius. This en- 
vironmental dependence of isophotal shape changes, if con- 
firmed, would imply that external environments have signif- 
icant effects on the inner properties of galaxies. 

There have been many theoretical interpretations for 
the origins of discy and boxy ellipticals. For example, Naab, 
Burkert & Hernquist (1999, also Nabb & Burkert 2003 and 
Naab, Khochfar & Burkert 2006) concluded that equal-mass 
mergers of two disk galaxies tend to produce boxy ellipticals 
while non-equal mergers produce discy ellipticals. Khochfar 
& Burkert (2005) used semi-analytical simulations to study 
the origin of the discy and boxy ellipticals. They concluded 
that the isophotal shapes of merger remnants depend not 
only on the mass ratio of the last major merger but also 
on the morphology of their progenitors and the subsequent 
gas infall. These theoretical studies and comparisons with 
high-redshift objects (Pasquali et al. 2006) all used the ob- 
servational constraints from the limited sample of Bender et 
al. (1989). Therefore, it is necessary to revisit the properties 
of isophotal shapes and their relations with other galactic 
properties in early-type galaxies using a larger sample and 
images with excellent qualities. 

The statistics of 04 /a parameter are important for a va- 
riety of applications, for example gravitational lensing. De- 
viations from the pure ellipses can induce higher-order sin- 
gularities such as butterflies and swallowtails, which can in 
turn produce sextuplet and octuplet images (Keeton, Mao 
& Witt 2000; Evans & Witt 2001). Evans & Witt (2001) 
estimated that perhaps ~ 1 per cent of galaxy-scale lenses 
are sextuplet and octuplet imaged systems. The study of 
a large and homogeneous sample of ellipticals and SO's is 
important for providing sufficient statistics of the shapes of 
E/S0 galaxies, which dominate the lensing cross-section on 
galaxy scales. 

In this paper, we use the Sloan Digital Sky Survey 
(SDSS) Data Release 4 (DR4) data of 847 nearby E/S0 
galaxies to study the distributions of their isophotal pa- 



rameters and the relations of isophotal shapes with other 
galactic properties. The SDSS survey provides high-quality 
photometry in five colours under typical seeings of about 1.5 
arcseconds over a large area of the sky. In §2, we discuss how 
we select our sample of elliptical and SO galaxies, outline the 
data reduction procedures and how we estimate the param- 
eters that describe the shapes of E/S0 galaxies. The size of 
our sample is about a factor of 8 larger than that used by 
Bender et al. (1989). This large sample size allows us to con- 
sider more carefully how the shapes of galaxies are correlated 
with a variety of parameters, such as their velocity disper- 
sions, luminosities, and colours. We also consider how their 
shapes are correlated with the environments. These results 
are presented in §3. We discuss our results in the last section 
(§4). Throughout this paper, we adopt a cosmology with a 
matter density parameter S7 m = 0.3, a cosmological constant 
JIa = 0.7 and a Hubble constant of Ho = 70kms -1 Mpc -1 . 



2 SAMPLE SELECTION AND DATA 
ANALYSIS 

2.1 Sample selection 

Our sample objects are selected from the SDSS DR4 pho- 
tometric catalogue. We select all objects according to the 
following criteria: (1) Their redshifts must be smaller than 
0.05; this is to ensure sufficient spatial resolutions to resolve 
these galaxies. (2) The velocity dispersions cover a range of 
200 kms -1 to 420 kms" 1 . The latter is approximately the 
maximum velocity dispersion that can be measured reliably. 
(3) The target galaxies should not be saturated or located 
at the edge of the corrected frame (Stoughton et al. 2002). 
In total, there are 1102 galaxies satisfying the above three 
criteria. However, 11 of these 1102 objects are unavailable 
from the SDSS DR4 Data Archive Server. Furthermore, to 
ensure the early-type galaxy morphology of our sample, we 
visually examined all images and find 916 (w 83 per cent) of 
them are E/S0 galaxies and not contaminated by companion 
galaxies or bright stars. To reduce the impact of dust on the 
photometry, objects with visible dust lanes are also excluded 
from our sample. As we will estimate the intensity- weighted 
Fourier coefficients and ellipticity, e, (see 12.31 within a re- 
gion from twice the seeing radius (2 r s ) to 1.5 times Petrosian 
half-light radius (1.5-Rso), objects with 2r s >= 1.5-Rso are 
excluded from the following analyses. The r-band seeing ra- 
dius for each object is the effective PSF width determined at 
the center of each frame by the SDSS pipeline (Stoughton et 
al. 2002). However, in most cases, our objects are not located 
at the centres. Stoughton et al. (2002) pointed out that the 
PSF FWHM can vary by up to 15% from one side of a CCD 
to the other, even in the absence of atmospheric inhomo- 
geneities. We tested the effect on isophotal parameters by 
enlarging the effective PSF width by 10% and found that 
the median differences for the weighted means of the 03/ 'a, 
a 4 /a, b z /a, b 4 /a and ellipticity are 0.072, 0.15, 0.020, 0.030 
and 2 times their corresponding errors respectively. The rel- 
atively large difference in ellipticity (~ 2a) between these 
two data sets is probably due to the monotonic ellipticity 
radial profiles for a large fraction of our sample objects and 
the small (formal) errors (a) in the ellipticity. 

For our sample objects, the redshifts are low and thus 
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should be corrected for the Local Group infall. We adopted 
the Local Group relative redshifts provided by Blanton et al. 
(2005) to calculate the luminosity distances. There are three 
objects without such information and another two with neg- 
ative redshifts. These five objects are excluded from the fol- 
lowing analyses. So our final catalogue includes 847 E/SOs. 



2.2 Data reduction 

2.2.1 Image reduction 

We use the SDSS r— band images in this paper to per- 
form isophotal analyses. The corrected frame obtained from 
SDSS is bias-subtracted, flat-fielded, and corrected for cos- 
mic rays and pixel defects (Stoughton et al. 2002). Before 
surface photometry, sky background fitting and subtraction 
are performed (Wu et al. 2005). To do this, we detect all ob- 
jects by SExtractor (Bertin & Arnouts 1996) and mask them 
to generate a background-only image. The sky background 
is obtained by fitting this image and then subtracted from 
the corrected frame. After background subtraction, the cor- 
rected frame is trimmed to a 501x501 frame which centres 
on the target. 

SExtractor was then run on the trimmed frame to gen- 
erate a "SEGMENTATION" image (See documents for SEx- 
tractor), which identifies all objects in the frame. A mask 
image with all detected objects except the galaxy of inter- 
est flagged can then be obtained from the "SEGMENTA- 
TION" image. We carefully examined all the mask images, 
and found that SExtractor performs satisfactorily for mask- 
ing neighbour objects in most cases. For a few (five) rare 
cases that SExtractor did not work well, we carefully tuned 
the parameters in SExtractor to create good mask images. 
Surface photometry is performed on the trimmed image us- 
ing the IRAF task ellipse (Jedrzejewski 1987) with the 
masked areas excluded from the fitting. This task requires 
some initial guesses for the geometric centre, ellipticity and 
position angle for the isophotes; these guesses are obtained 
by running SExtractor on the masked image. When run- 
ning ellipse on the targets, we allow the geometric centre, 
ellipticity and position angle to vary freely and assume a 
logarithmic step of 0.1 in radius along the semi-major axis. 



2.2.2 Isophotal shapes from Fourier analysis 

In the literature, the isophotes of elliptical galaxies have 
been measured in two different ways, which we outline be- 
low. 

(i) The first method first finds the isophote for a given 
surface brightness Jo (by interpolation, for example). The 
isophote can be described by R(9), where (R, 6) are the polar 
coordinates centred on the centre of an ellipse. In general, 
the isophote is not a perfect ellipse. The deviation of the 
isophote from a perfect ellipse can be expanded in Fourier 
series in the polar angle 

5R(6) = R(6) - R cl (6) = ao + ^2 

(a n cos n8 + b n sin n8) (1) 

where R c \(9) describes the best-fitting ellipse, and oo is the 
average deviation averaged over the polar angle. The first 
five coefficients (do, «i, 6i, 02 and 62) should be zero within 



the errors as they are found from fitting. The lowest order 
significant deviations are 03,63,04,64 etc. 

(ii) In the second method, an ellipse is drawn to approxi- 
mately match an isophote. The intensity along the ellipse is 
then expanded in Fourier series 

1(6) = I + ^2(A n cos n6 + B n sin n8) (2) 

where 7o is the intensity averaged over the ellipse, and 
A n and B n are the higher order Fourier coefficients. If 
an isophote is a perfect ellipse, then all the coefficients, 
(A n , B n ),n = 1, • • •, 00 will be exactly zero. 

The study by Bender et al. (1988) adopts the first approach 
while the IRAF task ellipse 1 adopts the second method. 
Ellipse outputs A n and B n divided by the semi-major axis 
length, a, and the local gradient. In the Appendix we show 
that this essentially gives a n /a and b n /a parameters, as used 
by Bender et al. (1988, 1989). So our results can be directly 
compared with those of Bender et al. (1988, 1989). 

2.3 Estimates of parameters 

The characteristic size (Petrosian half-light radius here), 
velocity dispersion and apparent Petrosian magnitude of 
galaxies are derived from the DR4 photometric catalogue. 
SDSS defines the Petrosian ratio IZp at a radius R from 
the centre of an object to be the ratio of the local surface 
brightness in an annulus at R to the mean surface brightness 
within R (Blanton et al. 2001; Stoughton et al. 2002): 

[ ( ^ R dR'2TTR'I(R')/[n(l.25 2 - 0.8 2 )R 2 } 

TIp(R) = v Jn v - — 3 

f* dR'2TvR'I(R')/(TvR 2 ) 

where I(R) is the azimuthally averaged surface brightness 
profile. The Petrosian radius _Rp is defined as the radius at 
which TZp(Rp) equals 0.2. The Petrosian flux in any band 
is then defined to be the flux within two Petrosian radii: 

f 2R P 

F P = / 2nR'dR'l(R'). (4) 
Jo 

The Petrosian half-light radius R50 is the radius within 
which 50 per cent of the Petrosian flux is contained. The 
Petrosian quantities provided by SDSS are derived using cir- 
cular apertures. For galaxies with large ellipticities, this will 
give a different R50 compared with that using elliptical aper- 
tures. However, the adoption of circular apertures will not 
affect our results statistically as there are few very elliptical 
galaxies in our sample (Fig.|3]l. To be more quantitative, we 
calculated the Petrosian R50 using elliptical apertures for 
some sample objects. Even for the most elliptical galaxies 
in our sample, the difference between circular and elliptical 
Petrosian R50 is only up to ~ 13%. Such a small change in 
R50 will not change our results significantly (see below). 

We derived the absolute Petrosian magnitude from the 
apparent Petrosian magnitude and the extinction (A) in 
each filter provided by SDSS DR4 photometric catalogue 
by M = m ~ 5 log(D L /10pc) - A. As the redshift of our 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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sample is low (< 0.05) and the range of the redshift is small 
(from 0.01 to 0.05), no k-correction is applied here. Using 
the k-corrections (corrected to z = 0.0) provided by Blan- 
ton et al. (2005; see also Blanton et al. 2003), we find that 
the median and maximum k-corrections in the r-band for 
our sample objects are 0.04 and 0.06 mag respectively, and 
hence this should not affect our results significantly. 

As the velocity dispersion provided by SDSS is that in 
the 3" fiber aperture, we correct it to a standard circular 
aperture defined to be one-eighth Petrosian half-light radius 
similar to Bernardi et al. (2003a). It should be noted that 
Bender et al. (1988, 1989) used the effective radius (half- 
light radius) from de Vaucouleurs law fitting as the char- 
acteristic size of a galaxy. We use the Petrosian half-light 
radius since it is a well defined quantity, while in contrast, 
the derivation of effective radius depends on how well the de 
Vaucouleurs law fits the surface brightness profile, the range 
of radius used in the fit (see Bender et al. 1989 and references 
therein) and how well the sky background is subtracted in 
the photometry. 

Isophotal surface photometry can provide the structural 
parameters, such as Fourier coefficients (a^/a, a^/a, bz/a, 
bi/a), ellipticity and position angle (PA) in each elliptic an- 
nulus, where a is the length of the semi-major axis of the 
fitted ellipse. 

In this paper, we mainly study isophotes using quanti- 
ties averaged over intensities or those at one or one and a 
half Petrosian half-light radii. Note that the radial profiles of 
isophotal parameters we used here are along the geometric 
mean radius, \fab~, where a and b are the lengths of semi- 
major and semi-minor axes of a fitted ellipse. To be more 
specific, the characteristic values of Fourier coefficients and 
ellipticity are calculated by weighting them with the inten- 
sity (counts) in the isophote and inversely with their RMS 
errors over a region of 2r s to 1.5-Rso- This is similar to that 
of Bender et al. (1988), who used the mean 0,3 /a and 0,4 /a 
between twice the seeing radius and 1.5 effective radii. Note, 
however, that Bender et al. (1988, 1989) used a different defi- 
nition of characteristic ellipticity - They adopted the peaked 
value for a peaked ellipticity radial profile and the ellipticity 
at the effective radius in case of a monotonically increasing 
or decreasing ellipticity radial profile. The changes in 03 /a, 
a,4,/a, ellipticity, PA and the centres of ellipses between R^o 
and 1.5i?5o are also estimated. The standard deviations of 
these quantities are calculated by error propagation. As ex- 
pected, fainter objects have larger errors than the brighter 
ones. All the data are available online 2 . 

The surface brightness profiles of E/SO's are often fitted 
with the Sersic profile (Sersic 1968). From the fitting, it is 
possible to obtain the sersic index n and a sersic effective 
radius within which half of the light is contained. Blanton 
et al. (2005) also gave the sersic index n and -R50 in their 
tables. We compared the SDSS Petrosian Rg,o and the sersic 
R50 given by Blanton et al. (2005) and found that they differ 
by 14 per cent on average. Blanton et al. (2005) argued that 
their sersic f?so is underestimated by ~ 10 per cent for large 
sizes and high sersic indices. We thus enlarged the Petrosian 
half-light radius Rg,o by 25 per cent and re-computed the 
isophotal parameters and compared them with those derived 
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by SDSS Petrosian f?so. The average differences between 
these two sets of data are 0.016, 0.023, 0.0045, 0.0073 and 
1.76 times the standard deviations for the weighted means of 
the a^/a, a^/a, bz/a, b^/a and ellipticity; the median differ- 
ence between these two data sets for the above five param- 
eters are respectively 0, 0, 0.0066a, 0.0049a and 0.67a, so 
our results are quite robust with respective to small uncer- 
tainties in .R50. The relatively large difference in ellipticity 
(average ~ 1.76a and median ~ 0.67a) is possibly due to 
the monotonic ellipticity radial profiles for a large fraction 
of our sample objects and the small (formal) errors (a) in 
ellipticity. 



3 RESULTS 

3.1 Sample properties 

Fig-0shows images for four of our sample galaxies. The top 
left and right panels show one nearly round and one very 
flattened ellipticals; both of which can be fitted by ellipses 
within noises. The bottom two panels show one example of 
a boxy elliptical and a discy elliptical respectively. It is clear 
that ellipticals and SO's can be reasonably fitted by ellipses 
but also show significant deviations. Below we first study 
some basic properties of our sample. 

The histogram of velocity dispersion (a) for our sample 
is shown in Fig. |5] The solid line shows the velocity disper- 
sion function derived by Sheth et al. (2003) from an analysis 
of a large sample of galaxies from the SDSS, carefully ac- 
counting for the selection effects and measurement errors. 
Our histogram matches the Sheth et al. curve reasonably 
well except at the small dispersion end, where our sample 
seems to show a small deficit but still within error bars. 
The reasonable match between our sample and the Sheth 
et al. curve indicates our sample does not have a strong 
bias regarding to the velocity dispersion. As the luminos- 
ity is correlated with the velocity dispersion through the 
Faber- Jackson relation (Faber & Jackson 1976; Bernardi et 
al. 2003b), we do not expect our sample to be significantly 
biased with respect to the luminosity (or magnitude). 

The ellipticity distribution for our sample is shown in 
Fig. [3] The distribution shows a peak around e « 0.2 and 
drops off sharply above e = 0.5 (correspond to an E5 ellipti- 
cal). The mean and dispersion in the ellipticity are 0.23 and 
0.13 respectively. Our distribution appears to show smaller 
ellipticities than the one shown in Lambas, Maddox & Love- 
day (1992) based on the APM survey which uses photo- 
graphic plates. Their distribution peaks around e = 0.25 
and drops off sharply above e = 0.8. As they adopted a dif- 
ferent definition of ellipticity (their equations 1 and 2), we 
used similar quantities as provided by SDSS (See Kuehn & 
Ryden 2005 for details) to examine the ellipticity distribu- 
tion. The distributions of these two ellipticities in the SDSS 
data are very similar. Therefore, the difference between our 
sample and theirs is not due to different definitions of ellip- 
ticity; it is not clear what causes it. However, our ellipticity 
distribution is in very good agreement with Fasano & Vio 
(1991), who used CCD images of elliptical galaxies. The dif- 
ferential probability distribution of the ellipticity can be well 
fitted by a Gaussian: 
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Figure 1. Four examples of our sample galaxies. The top two panels show two ellipticals that are well fitted by ellipses within noises. 
The ellipticities of these two galaxies are respectively 0.03 and 0.46. The bottom two panels show one boxy (left, cn/a = —0.02) and one 
discy (right, a^/a = 0.035) elliptical galaxy. The field size for each image is 34.4" by 34.4". Five isophotes are shown together with the 
best-fitting ellipses. 
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where the dispersion <r e = 0.14. This fit is shown as the 
solid curve in Fig. |31 The ellipticity distribution (dashed 
histogram) of the lens galaxy sample of Koopmans et al. 
(2006) is shown for comparison in the same Figure. The 
Kolmogorov-Smirnov test indicates a probability of 58.6 per 
cent that the two samples are drawn from the same parent 
distribution (for more see the discussion). 

The histograms of a^/a, a±/a, b^/a and 64/a are shown 
in Fig. [1] The distributions of 03/0, bz/a and 64/a are rea- 
sonably fitted by a single standard Gaussian with zero means 
and dispersions of 0.0032, 0.0035 and 0.0039 respectively. 
The histogram for 04/05, on the other hand, shows more 
asymmetries, particularly toward the positive a^ja indicat- 
ing an excess of discy E/SOs compared to boxy E/SOs. The 
distribution of 04/0 can be reasonably fitted by the sum 
of two Gaussian functions (The normalized \ 2 P er degree 



of freedom is 0.81, where the errors are taken as Poisson 
errors) : 

dP 4.61 x TO" 4 -A 4.32 x TO" 4 - <^-°-°°™) 2 

— - = = e 2 "i + — e 2 "2 (6) 

aAn \/2irai \/2ttct2 

where Aa = 04/0, o\ = 0.0035, and 02 = 0.01. 



3.2 Correlations with galaxy properties 

Fig. |S] shows how a^/a varies as a function of the mean 
ellipticity, e. It is clear that when e > 0.4, the isophotes show 
significant systematic deviations from perfect ellipses, with 
the a^ja parameter reaching 0.08 or so. As the ellipticity is 
both a function of the intrinsic shape and the viewing angle 
along the line of sight, this trend indicates that E/SOs, when 
viewed edge on, are more likely to show significant deviations 
from perfect ellipses. 

Fig. |S| shows how the as /a and a^ja parameters vary 
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Figure 2. Histogram of velocity dispersion without aperture correction for our 847 galaxies. The solid line is the velocity dispersion 
function derived by Sheth et al. (2003) from the SDSS. The error bars indicate the Poisson errors. For all histograms in this paper, the 
error bars have the same meaning and this will not be stated in the following. 



as a function of the velocity dispersion, a. There is no sig- 
nificant trend for the 03 /a parameter. For the 04 /a param- 
eter, there seems to be a slight indication that at small ve- 
locity dispersions (a < 250 kms" 1 ), there are more E/SOs 
with significant non-zero 0,4/0, parameters, around 0.05. This 
indicates a weak trend that more massive ellipticals may 
have less significant deviations from perfect ellipses. The 
trend of 03/0 parameter as a function of r-band absolute 
magnitude is similar to that of 03 /a versus a. The cor- 
relation between 04/0 and the r— band absolute magni- 
tude is more significant (Fig. |7J . The Spearman- Rank Or- 
der correlation coefficient and the probability that no cor- 
relation exists for 03/0 vs. a, 04/0 vs. a, 03/0 vs. M r , 
and 04/0 vs. M r are r s = -0.019,-0.042,0.0054,0.32 and 
Prob = 0.58,0.22,0.88,1.84 x 10" 21 respectively. Among 
these correlations, the relation between 04/0 and M r , the 
r— band absolute magnitude, is most significant; the frac- 
tion of boxy galaxies increases with increasing luminosity. 
The relations between 03/0, 04/0 and the dynamical mass, 
defined as 9a 2 R 50 /G (Bender et al. 1989), are similar to 
those with the r— band absolute magnitude. 

Faber et al. (1997) concluded that galaxies with My < 
—22 have 'core' luminosity profile, those with My > —20.5 
have 'power-law' luminosity profile and at intermediate 
(with —22 < My < —20.5) magnitudes, 'core' and 'power- 
law' coexist. Many studies (e.g., Rest et al. 2001) indicated 
that 'core' galaxies are more luminous, slowly rotating and 
boxy systems while 'power-law' galaxies are less luminous, 
fast rotating and discy systems. Therefore, it is interesting 



to estimate quantitatively the relative number of discy and 
boxy E/SOs in these magnitude intervals. We use the pre- 
scription of V = g — 0.55(g — r) — 0.03 given by Smith et al. 
(2002) to transform the (?-band magnitude into the l/-band 
magnitude for our sample objects. Using our isophotal anal- 
ysis results, we find that the ratios of the number of discy 
to boxy E/SOs in these magnitude intervals (My > —20.5, 
-22 < M v < -20.5 and M v < -22) are 2.8 (28/10), 
2.09 (466/223) and 1.07 (62/58) respectively. These num- 
bers show that the largest fraction of boxy E/SOs is in the 
brightest subsample, while the number of discy E/SOs in- 
creases relative to boxy E/SOs at the faint end. 

Fig.|H]shows the mean ellipticity, e, as a function of the 
r— band absolute magnitude and velocity dispersion. The 
data points with error bars show the median values and 
the lower (25 per cent) and upper (75 per cent) quartiles 
of e in bins of absolute magnitude (left panel) and veloc- 
ity dispersion (right panel). In many places in this paper, 
we quote the median and lower and upper quartiles as they 
are less sensitive to outliers (we sometimes also quote the 
usual mean and standard deviation value). As can be seen 
from the left panel, the ellipticities of E/SOs are correlated 
with the absolute magnitudes, although with large scatters: 
E/SOs become rounder with increasing luminosities, which is 
consistent with Vincent & Ryden (2005). In contrast, the el- 
lipticities of E/SOs show no significant correlation with their 
velocity dispersions. Quantitatively, a Spearman-Rank Or- 
der correlation analysis gives a correlation coefficient of 0.26 
at > 99.9 per cent significance level for e versus r— band 
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Figure 3. Histogram of ellipticity, e, for our sample (solid curve), e is the average ellipticity weighted by intensity between twice the 
seeing and 1.5 Petrosian half-light radii. The solid line shows an empirical Gaussian fit to the data. The dashed histogram shows the 
ellipticity histogram for the 15 lenses studied by Koopmans et al. (2006). 



absolute magnitude, while for e versus velocity dispersion, 
the correlation coefficient is 0.035 at a significance level of 
69.2 per cent. We also find that the mean ellipticity has a 
strong correlation with the dynamical mass, with a corre- 
lation coefficient of —0.27 at > 99.9 per cent significance 
level. 

Fig. |5] shows as/a and a^ja versus the g — r colour. 
Neither shows any strong correlation. Furthermore, Fig. 1101 
shows that the mean ellipticities of galaxies also have virtu- 
ally no correlation with their colours. 

It appears somewhat puzzling why both the mean el- 
lipticity and 04 1 a are significantly correlated with r-band 
absolute magnitude while they do not correlate with the 
colour or velocity dispersion, as if ellipticals/SO's follow the 
colour-magnitude relation and the Faber- Jackson relation 
(L oc a 4 ). However, these may be explained by the large 
scatters in these correlations - mean ellipticity versus M T 
(Fig. |HJ, a^ja versus M T (Fig. |7J, colour- magnitude rela- 
tion (e.g., Chang et al. 2006) and the Faber- Jackson relation 
(e.g., Bernardi et al. 2003b). Our sample selection criterion 
may also be partly responsible for the absence of correla- 
tions of the mean ellipticity and a^/a with velocity disper- 
sion since our velocity dispersion is restricted to within a 
factor of two ( H2.H . 

Ryden et al. (2001) found that the ages of elliptical 
galaxies are correlated with their ellipticities and luminosity 
profiles ('core' or 'power-law'). We cross-correlate our sam- 
ple objects with DR4 catalogue generated by Kauffmann et 
al. (2003) and find 345 objects in common. They derived 



the age using a Bayes method by exploring the D n (4000), 
H<5a and a library of star formation histories generated by 
the stellar population synthesis code provided by Bruzual & 
Chariot (2003; See Kauffmann et al. 2003 for details). We 
used here the median estimates of the r-band weighted mean 
stellar age in their table. However, no strong correlations are 
found between the age versus ellipticity and a^/a. 

As mentioned before, the sersic index n from surface 
brightness fitting is available from Blanton et al. (2005). We 
correlated this parameter with other isophotal parameters. 
No correlation was found between n and as/ a. However, 
there appears to be a correlation between n and a^/a: the 
Spearman-Rank Order correlation coefficient is —0.104 and 
the probability that no correlation exists between these two 
parameters is only 0.24 per cent. Galaxies with large sersic 
index n tend to have smaller deviations from perfect ellipses. 
There is also a weak trend that more flattened (i.e. larger 
ellipticity) galaxies have smaller values of n, but the scatter 
is quite large. 

3.3 Isophotal twists 

It is well known that isophotes in elliptical galaxies are not 
concentric ellipses as a function of radius. In this subsection, 
we concentrate on the changes of isophotes from one to one 
and a half Petrosian half-light radii (R50 — * 1.5i?5o). 

Fig. 1111 shows the change of ellipticity as a function of 
the velocity dispersion, colour, mean ellipticity, and the dif- 
ference in the position angle between R$o and 1.5 R$o- There 
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Figure 4. Histograms of 03 /a (upper left), (14/12 (upper right), bz/a (bottom left), and 64/12 (bottom right). The solid lines are the 
best-fitting single Gaussian model for 03 /a, 63 /a, 64/11 and double Gaussian model for (24 /a (See text). 



is no strong correlation between the change and either veloc- 
ity dispersion, colour or the change in position angle. There 
is however a correlation between the mean ellipticity e and 
the ellipticity change, with Spearman-Rank Order correla- 
tion coefficient r s = —0.373 at a significance level of > 99.9 
per cent. Elliptical galaxies with higher ellipticities tend to 
have a more negative value of Se — e(Rso) — e(1.5i?5o). This 
implies that the isophotes in the most elliptical galaxies be- 
come more flattened as one goes to larger radii. 

Fie, ll2l shows the histogram of the ellipticity change, Se. 
The mean change in ellipticity is about —0.03 and the scatter 
is about 0.051. The differential probability distribution can 
be well-fitted by the sum of two Gaussians: 

dP 0.017 - is '- 5 ? )2 0.002 - <s "- li '? )2 

= t= e 2 "i + e 2 '2 (7) 

dOe V27TCT1 V27TCT2 

where <5ei = -0.018, o\ = 0.035, Se 2 = -0.124, and a 2 = 
0.034. This fit is shown as the solid curve. The second peak 
is due to the galaxies with high ellipticities (see the bottom 
left panel of Fig. ITU. 

Fig. 1131 shows the change of position angle as a func- 
tion of the velocity dispersion, colour, mean ellipticity, and 
the difference in ellipticity between R50 and 1.5-Rso. In this 
case, there is no strong correlation with all these quantities. 
Fig. 1141 shows the histogram of the changes of the position 
angles. The differential probability distribution can be well 
fitted by a power-law 

dP , a ,,/", 5PA\~ 13 



where <5PA = 1.07° and /3 = 1.81. The median value of 
SPA is 1.61° . In contrast, the mean change in position an- 
gle is much larger, about 4.12° and the standard deviation 
is about 8.39°. The large standard deviation is due to an 
extended tail in SPA. 

Figs. 1151 and [Hal show the changes of the centres of el- 
lipses along the major (x) and minor (y) axes as a function 
of velocity dispersion, colour, mean ellipticity and the dif- 
ference in position angle between R50 and 1.5iZso- There 
is no significant correlation among these quantities. The 
histograms of the changes of the centres of ellipses are 
shown in Fig. 1171 The differential probability distribution 
can be fitted by a standard Gaussian with zero means and 
c«x/(0.5flB ) = 0.030 and a Sy/(0 , 5 Rso) = 0.021 respectively. 
The centroid shifts are hence small, only a few percent of 
the Petrosian half-light radius. 

3.4 Radio properties of elliptical galaxies 

We search the radio counterparts of our early-type galaxies 
using the FIRST (Becker, White & Helfand 1995) and NVSS 
(Condon et al. 1998) catalogues. The NVSS and FIRST sur- 
veys provide an angular resolution (beam size) of 45" and 
5" respectively. However, the centroids can be determined 
much more accurately and the accuracy scales roughly as 
1/V s ig ria l — to — noise ratio. We adopt a conservative cross- 
correlation radius between the radio catalogue and our sam- 
ple objects. We found 162 counterparts for FIRST within a 
5" cross-correlation radius and 130 for NVSS galaxies within 
a 20" cross- correlation radius. 
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Figure 5. 04/0 as a function of cllipticity. The data points with error bars are the median, lower (25 per cent) and upper (75 per cent) 
quartiles for galaxies in bins of width 0.1 in ellipticity except the last bin with a width of 0.3 to include all the remaining objects. The 
solid curve shows the boundary between quadruple and sextuple lenses derived by Evans & Witt (2001). The median error bars are 
shown at the top right corner. All error bars shown at the top right corner of the following figures have the same meaning. 



As the position accuracy for FIRST is higher, here we 
only present the results for the objects cross-correlated from 
FIRST and SDSS. Fig. [TH] shows the 1.4GHz FIRST peak 
luminosity as a function of the as /a parameter and 0,4 /a 
parameter. There is no trend between the radio power and 
the 0,3 /a parameter. For the 0,4/ a parameter (shown in the 
right panel), there appears to be a lack of galaxies at the top 
right corner. (The trend is similar for the 1.4GHz FIRST in- 
tegrated luminosity versus a^/a.) Thus we confirm the trend 
found by Bender et al. (1989) that discy ellipticals (with 
0,4 > 0) are weaker radio emitters than their boxy counter- 
parts, although there are much less fainter radio emitters in 
our sample compared to Bender et al. (1989) because of the 
larger redshift in our sample relative to theirs. 

3.5 Correlations with environment 

Di Tullio (1979) found that the ellipticity of isolated ellip- 
ticals on average decreases as a function of radius, while 
for ellipticals in clusters and groups, the ellipticity can 
either increases, decreases or peaks. This conclusion was 
based on photographic plates for 75 ellipticals in clusters, 
groups, pairs and isolated galaxies. More recently, using the 
SDSS data, Kuehn & Ryden (2005) concluded that as the 
local density of galaxies increases, galaxies with de Vau- 
couleurs profiles become rounder in the magnitude range 
of — 22 < M r < —18. On the other hand, for the most lumi- 
nous galaxies (M r < —22), galaxies become rounder in the 



inner regions, but exhibit no change in their outer regions as 
the local density of galaxies increases. Shioya & Taniguchi 
(1993) used 96 galaxies in Bender et al. (1989) to study the 
isophotal shape-environment relation and found that boxy- 
type galaxies favor the dense environments while discy-type 
galaxies prefer the field environments. 

In this paper, following Kuehn & Ryden (2005), we use 
a measure of the local density as the number of galaxies 
within a cylinder in the redshift space. The depth of the 
cylinder (centred on a target) is 12/i _1 Mpc, and we vary the 
radius of the cylinder from 0.3/i -1 Mpc to 0.5/i -1 Mpc, where 
h = #o/(100 kms -1 Mpc) = 0.7. According to Kuehn & Ry- 
den (2005), if a target galaxy is within one cylinder width 
of the SDSS spectroscopic survey border, it is removed from 
the analyses. The radii of the cylinders we adopted above 
are smaller than those used by Kuehn & Ryden (2005) be- 
cause if we adopt their larger values (e.g., 2h~ Mpc), too 
few galaxies would be left because our sample galaxies are 
at lower redshift and so the boundary effects discussed above 
would exclude too many galaxies. Following Kuehn & Ryden 
(2005), we define a volume limited sample which is selected 
from SDSS DR4 catalogue by constraining the spectral red- 
shift smaller than 0.055. For these objects, we also adopt 
redshifts corrected for the Local Group infall provided by 
Blanton et al. (2005). As the maximum redshift of our sam- 
ple galaxies is 0.05 and the limiting apparent magnitude 
of the spectroscopic survey is m r = 17.77, the correspond- 
ing absolute magnitude is M t = —19. We use the number 
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Figure 6. as /a (left) and m/a (right) as a function of velocity dispersion. The data points with error bars are the median, lower (25 
per cent) and upper (75 per cent) quartiles for 5 bins separated by 20kms — 1 from <r = 200, 300 km s" 1 . The last bin has a large bin 
width of 50kms _1 to include all the remaining galaxies. 



counts of galaxies with M T < — 19 within a cylinder of the 
target galaxy as an indicator of the local density. This is 
only an approximate measure of the local environment, due 
to the effects of fiber collisions (see Kuehn & Ryden 2005 for 
detailed discussions) and redshift distortions. But it is suffi- 
cient for our purpose. For brevity, we term galaxies without 
neighbours in the cylinder as 'isolated' galaxies. 

Fig. 1191 shows the 04/0 distribution for 'isolated' galax- 
ies (solid curve) and galaxies with neighbour (dot-dashed 
curve). The left and right panels are for a cylinder with 
radius of 0.3/i _1 Mpc and 0.5/i _1 Mpc respectively. We can 
see that the fraction of discy galaxies is larger in isolated 
environments. Quantitatively, for cylinders with radii of 
0.3/i" 1 Mpc (left panel) and 0.57l _1 Mpc (right panel), the 
ratios of the number of discy to boxy E/SOs in 'isolated' en- 
vironments are 3.14 (138/44) and 2.92 (76/26) respectively. 
In contrast, for galaxies with neighbours, the correspond- 
ing ratios are 1.63 (381/234) and 1.72 (386/225). Clearly 
discy E/SOs are more abundant in isolated environments, 
in agreement with the conclusion of Shioya & Taniguchi 
(1993). The Kolmogorov-Smirnov test indicates a probabil- 
ity of 99.9 (90.2) per cent that the two samples are different 
for cylinders with radii of 0.3/i _1 Mpc (0.5/i _1 Mpc). 

The change of ellipticity as a function of the local den- 
sity is shown in Fig. 1201 There seems no correlation between 
the change of ellipticity and the local density. This result is 
not consistent with Di Tullio (1979). His sample based on 
photographic plates may be too small to be statistically con- 
vincing. Reda et al. (2004) presented the ellipticity profiles 



for 14 isolated early- type galaxies (their fig. 4); their re- 
sults also do not support the conclusion of Di Tullio (1979). 
Fig. 1211 shows the mean ellipticity as a function of the local 
density, a plot similar to figs. 6 and 7 in Kuehn & Ryden 
(2005). Qualitatively, our results are consistent with Kuehn 
& Ryden (2005). 

As isophotal twists may be caused by external pertur- 
bations, we examined the isophotal twist - local density re- 
lation in Fig. 1221 However, there is no correlation found. 
Alternatively, we also examined the ratio of the number of 
E/SOs with large PA changes between R50 and 1.5i?so to 
that of E/SOs with small PA changes for isolated galaxies 
and galaxies with neighbours. The results show that the 
fraction of E/SOs with large PA changes is slightly larger 
for galaxies with neighbours than that for isolated galaxies. 



4 CONCLUSIONS AND DISCUSSIONS 

In this paper, we studied the shapes of isophotes for a sam- 
ple of 847 nearby E/S0 galaxies selected from the SDSS 
DR4 catalogue. We derived the distributions of the intensity- 
weighted ellipticity and Fourier coefficients and also studied 
the changes in the isophote shapes between one and one 
and a half Petrosian half-light radii. For many of these dis- 
tributions, we provide empirical fitting formulas of the data, 
which may be useful for comparisons with theoretical stud- 
ies. We further examined many possible correlations between 
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Figure 7. as /a (left) and a^/a (right) as a function of the absolute magnitude in the r-band. The median, lower (25 per cent) and upper 
(75 per cent) quartiles for all galaxies are shown in 5 bins with bin width of 0.5 mag except the first and last bin which have a width of 
1.5 and 1 mag respectively. 



the isophote shapes and physical properties of E/SO galax- 
ies. Our main conclusions are the following: 

(i) The mean ellipticity of galaxies is about 0.23 with a 
scatter of 0.13. This is somewhat narrower than the distribu- 
tion found by Lambas et al. (1992). The intensity- weighted 
mean ellipticity shows no strong correlation with either the 
velocity dispersion or the g — r colour. However, it is corre- 
lated with the r— band absolute magnitude and dynamical 
mass. These correlations indicate that the more luminous 
(massive) E/SOs are rounder than the fainter (less massive) 
ones. 

(ii) The most significant deviation of the isophotes from 
ellipses is the a^ja parameter determined from Fourier se- 
ries. This parameter shows no significant correlation with 
colour and velocity dispersion. It is, however, correlated with 
the r— band absolute magnitude, in agreement with Bender 
et al. (1989). As the luminosity decreases, the number of 
discy E/SO's increases relative to boxy ones. The a^/a pa- 
rameter is also strongly correlated with the mean elliptic- 
ity (Fig. E/SO's with higher ellipticities tend to show 
stronger deviations in isophotes from perfect ellipses. They 
can become either more discy or more boxy. This suggests 
that the isophotes are likely a function of viewing angle 
- edge-on ellipticals will show stronger deviations from el- 
lipses. 

(iii) Our studies on the relation between isophote shapes 
and environments show that discy E/SOs favor field envi- 
ronments while boxy E/SOs prefer dense environments, in 



agreement with Shioya & Taniguchi (1993). The changes of 
ellipticity and position angle with radius, however, have no 
significant correlations with the environments. 

(iv) We confirm the results of Bender et al. (1989) 
that discy ellipticals tend to lack powerful radio emitters 
(Fig. I18H . although the trend is not dramatic. 

(v) The changes in ellipticity and position angle between 
one and one and a half Petrosian half-light radii show no 
correlations with galaxy properties, except that the change 
in ellipticity is correlated with the mean ellipticity (see 
Fig. llH , E/SO galaxies with larger ellipticities become more 
flattened at large radii. The median change in the ellipticity 
and position angle is about —0.023 and 1.61° respectively. 
The change in the position angle has a significant tail. 

Gravitational lensing can be used to probe the shape of 
total matter distribution in intermediate redshift galaxies. 
Yoo et al. (2006) derived the limits for four lensing sys- 
tems, SDSS J0924+0219, HE 0435-1223, B 1938+666, and 
PG 1115+080. The upper limits on a 3 /a are 0.023, 0.019, 
0.037 and 0.035, and a 4 /a are 0.034, 0.041, 0.051, 0.064 
respectively for these galaxies. Our study can only con- 
strain the shape of stellar matter distribution. Neverthe- 
less, Fig. 01 shows that the vast majority of our elliptical 
galaxies satisfy these constraints as long as the dark matter 
distribution does not show more significant deviations from 
the elliptical distributions. It is interesting to note that in 
Fig- El the ellipticity distribution (dashed histogram) of the 
lens galaxies (selected largely by mass) are consistent with 
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Figure 8. Ellipticity vs. r-band absolute magnitude (left) and velocity dispersion (right) for our sample galaxies. In the left (right) 
panel, the data points with error bars have the same meaning with Fig. 171 (Fig. 
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Figure 10. Mean ellipticity, e, as a function of color, g — r. 



our luminosity-selected sample. The number of lenses in the 
Koopmans et al. study is only 15, so is still somewhat small 
to draw stronger conclusions. Another effect of gravitational 
lensing is that deviations from perfect ellipses can cause un- 
common lensing configurations, such as sextuplet and oc- 
tuplet imaging. The solid curves in Fig. Q>] show the critical 
lines for forming sextuplet lenses (eq. [32] in Evans & Witt 
2001). There are only 11 galaxies that fall outside the criti- 
cal lines (and one almost on the line). Hence if mass roughly 
follows light and there is no significant evolution in shapes 
of ellipticals as a function of redshift (see Pasquali et al. 
2006), then the fraction of galaxies that can produce more 
than 4 images will be of the order of ~ 1 per cent. However, 
the magnification bias (e.g., Turner, Ostriker & Gott 1984) 
for these systems is likely somewhat higher, bringing the 
fraction somewhat higher. So far, about 100 gravitational 
galaxy-scale lenses are known, and only one system shows 
six images due to a compact group of galaxies (B 1359+ 154, 
Rusin et al. 2001). Currently the number of known lenses 
(~ 100) is too small to allow a definitive comparison. 

Hydrodynamical simulations of galaxies cannot yet sim- 
ulate individual elliptical galaxy realistically (see, e.g. Meza 
et al. 2003 for an attempt), let alone the population as a 
whole. The statistical results we obtained in this paper can 
serve as a comparison sample for future numerical simula- 
tions. An intermediate, and profitable step is to use collision- 
less simulations of disk galaxy mergers to study the isopho- 
tal shapes (Naab & Burkert 2003; Jesseit, Naab & Burkert 
2005). Naab & Burkert (2003) ran 112 merger simulations 
with mass ratios of 1:1, 2:1, 3:1 and 4:1. Their 1:1 disk galaxy 
mergers can reproduce the boxy ellipticals, but most of the 



simulated ellipticals have cu/a < 0.03 whatever the mass 
ratio (see their fig. 6), so it will be difficult to reproduce the 
most discy E/SOs {a^/a >= 0.03), about 5.9 per cent of our 
sample galaxies. Their origin remains unclear. Perhaps the 
inclusion of dissipation can remedy this situation. 

Weil & Hernquist (1996) argued that pair and multi- 
ple merger remnants (e.g., from groups) show differences in 
both spatial and kinematic properties. In particular, they 
found that multiple merger remnants are more likely to ap- 
pear nearly round from many viewing angles. If the multiple 
mergers produce more luminous ellipticals, then this would 
predict a correlation between the luminosity and the ellip- 
ticity. Fig.[S]shows that there is indeed some correlation be- 
tween the luminosity and dynamical mass versus ellipticity. 
It would be interesting to further explore these correlations 
with more simulations. 
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Figure 12. Histogram of the difference in the ellipticity between one and one and a half Petrosian half-light radii for our sample galaxies. 
The solid line is a double Gaussian fit to the data. 
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Figure 13. The difference in the position angic between one and one and a half Petrosian half-light radii for our sample galaxies vs. 
velocity dispersion, colour, mean ellipticity, and change in the ellipticity between R50 and 1.5 Rso- 



APPENDIX A: FITTING ISOPHOTES WITH 
ELLIPSES 

Let us consider the surface brightness on the best-fitting 
ellipse R e i(9) given in Eq. Q. If 8R is small, then we have 
for a fixed 8 



I(Rei(0)) = I(R - SR) « I(R) 
Therefore, 



dl_ 

dR 



5R 



I(R) = I(R cl {9)) + 



dl_ 

dR 



6R 



(Al) 



(A2) 



As (I(R el (6))) w / , we have 

[ a ° ( a ™ cos n ^ ^™ s * n (A3) 



«R)=i + § i 



where we have used the Taylor expansion to the first order 
in the second step. By definition, the right hand side of the 
equation must be equal to that of Eq. @. From this, we 
obtain 



dl_ 

dR 



The above equation implies 

a 70 



(A4) 



(A5) 



where a is the semi-major axis length and we have used 



dl_ 

dR 



(A6) 



is the local radial intensity gradient. The ellipse program 
in IRAF outputs the right hand side of Eq. 1A5L which is 
identical to the parameter used in Bender et al. (1989) - the 
left hand side of Eq. IA5I . 



Isophotal Shapes of E/SO galaxies 17 




10 20 30 



|PA(R 60 )-PA(1.5RJ| (deg) 

Figure 14. Histogram of the difference in position angle between one and one and a half Petrosian half-light radii for our sample galaxies. 
The solid line is a power-law fit to the data; the power-law index is —1.81. 
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Figure 15. Controid shift (in units of half a Pctrosian half-light radius) from one to one and a half Petrosian half-light radii in the 
major-axis direction vs. velocity dispersion, colour, mean ellipticity, and change in the position angle between R50 and 1.5i?5o- 
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Figure 16. Centroid shift (in units of half a Pctrosian half-light radius) from one to one and a half Petrosian half-light radii in the 
minor-axis direction vs. velocity dispersion, colour, mean ellipticity, and change in the position angle between R50 and 1.5 -Rso- 
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Figure 17. Histograms of the centroid shifts from one to one and a half Petrosian half-light radii in the major (left) and minor (right) 
axis directions in units of half a Petrosian half-light radius. The solid lines are the best-fitting Gaussians. 
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Figure 18. Radio power vs. a-.il a and a^/a parameters for elliptical and SO galaxies identified from a cross-correlation between the SDSS 
and the FIRST radio catalog. 
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Figure 19. Histograms of a^/a for isolated galaxies (solid curve) and galaxies with neighbour (dot-dashed curve). Isolated galaxies are 
those with zero neighbour within a cylinder of depth 12/i -1 Mpc and radius 0.3/i _1 Mpc (left) and 0.5/i — 1 Mpc (right). For clarity, Poisson 
errors are not plotted. 
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Figure 20. The change of ellipticity between one and one and a half Petrosian half-light radii as a function of the local density, as 
measured by the number of galaxies, AT cy i inc j er , within a cylinder of depth 12h -1 Mpc and radius 0.3/i -1 Mpc (left) and 0.5ft _1 Mpc 
(right). The data points with error bars are the median, lower (25 per cent) and upper (75 per cent) quartiles for bins of AT cy i inc j er . 
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Figure 21. Ellipticity versus local density, as measured by the number of galaxies within a cylinder of depth 12/i _1 Mpc and radius 
0.3/i _1 Mpc (left) and 0.5/i _1 Mpc (right). The data points with error bars are the median, lower (25 per cent) and upper (75 per cent) 
quartiles for bins of N C y\i n d CI ■ 
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Figure 22. The twist of position angle from one to one and a half Petrosian half-light radii as a function of local density, as measured 
by the number of galaxies within a cylinder of depth 12/i _1 Mpc and radius 0.3h -1 Mpc (left) and 0.5/i _1 Mpc (right). The data points 
with error bars are the median, lower (25 per cent) and upper (75 per cent) quartiles for bins of AT cylillder . 



